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Abstract 

We study one-dimensional models of particle diffusion and attachment /detachment from islands 
where the detachment rates 7 of particles at the cluster edges increase with cluster mass m. 
They are expected to mimic the effects of lattice mismatch with the substrate and/or long-range 
repulsive interactions that work against the formation of long islands. Short-range attraction is 
represented by an overall factor e <C 1 in the detachment rates (masses m > 2) relatively to 
isolated particle hopping rates (e ~ e~ E l T , where E is a binding energy and T is the temperature). 
We consider various mass-dependences of 7(771), from rapidly increasing forms such as 7(771) ~ m 
to slowly increasing ones, such as 7(771) ~ , with constant b. A mapping onto a column 

problem shows that these systems are zero-range processes, whose steady states properties are 
exactly calculated under the assumption of independent column heights in the Master equation. 
Simulation provides accurate island size distributions which confirm analytic reductions and are 
particularly useful whenever the analytical tools cannot provide results in closed form. The shape 
of island size distributions can be changed from monomodal to monotonically decreasing by tuning 
the temperature or changing the coverage (one-dimensional density p). In all cases, small values 
of the scaling variable X = e~ 1 p/ (I — p) favour the monotonically decreasing ones. However, 
for large X, rapidly increasing rates 7(771) lead to distributions with peaks very close to (m) and 
rapidly decreasing tails, while slowly increasing 7(777) provide peaks close to (m)/2 and fat right 
tails. 
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I. INTRODUCTION 



In order to understand the main processes that take place during the growth of thin 
films and multilayers, it is important to study the initial steps of those processes, i. e. 
the submonolayer regime. The current belief is that insight derived from the combination 
of experiment and modeling can lead to better control of nanostructures formed during 
deposition. Consequently, models of submonolayer growth and coarsening with several types 
of interactions between the adparticles were intensively studied in the last decades - for a 
recent review with applications to homoepitaxial growth, see Ref. jl|. Non-equilibrium 
statistical models are particularly useful for modeling real systems that are confined to 
long-life metastable states while some atomic processes take place in next-to-equilibrium 
conditions with the environment. 

The recent advances on the production of elongated structures (e. g. nanowires) along 
st ep edges of vicinaa surfaces Jfl fl fl Q and other highly a» ls o tr opie growth processes 
motivated the study of models thatproduce effectively one-dimensional structures on sur- 



faces under certain conditions 



s tnatproc 



11] . However, in many cases it is advantageous to 



restrict the processes of particle diffusion and attachment/detachment from islands (possibly 
competing with the continuous atom deposition) to one spatial dimension, since it enables 
a more detailed study of the effects of different physico-chemical parameters. One of the 
interesting problems in this field is to determine the island size distributions under different 
growth conditions. For instance, recent works have debated the mechanisms responsible for 
the onset of monomodal or monotonically decreasing distributions in submonolayer growth 
on vicinal surfaces, since both types of distribution were already observed experimentally 

The peaked island size distributions are the most frequently observed ingrowth on planar 



surfaces [l| and also appear in some effectively one-dimensional systems (3, E], |l2j] . In theo- 
retical models, this feature is observed when some critical island size is assumed, above which 
there is no further atom detachment. This occurs both for deposition competing with diffu- 
sion and for post-deposition coarsening. On the other hand, statistical equilibrium models, 
where reversibility of attachment and detachment from islands is implicitly assumed, provide 
monotonically decreasing size distributions in one dimension (more precisely, exponentially 



decreasing ones) 



0,Q 
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The simplest reversible non-equilibrium models consider the processes illustrated in Fig. 
a in its one- dimensional version (after deposition has stopped, i. e. with conserved mass) 
"bl [l5 . Q, Q. They account for particle diffusion, aggregation and detachment from 
islands, with no critical island size. For simplicity, it is assumed that attachment occurs 
immediately upon contact of a diffusing particle with a cluster, while the detachment occurs 
with rate 7(771), where m is the mass of the cluster from which the particle leaves. Previous 
work in one-dimension only considered the case of constant 7 for all islands (m < 2) 19], 
where the steady state presents monotonically decreasing cluster size distributions, similarly 
to the equilibrium models 0, [3]. The great applicability of the two-dimensional versions of 
these models explains the small number of studies of the one-dimensional cases. 

In this paper, we will consider this class of non-equilibrium one-dimensional models with 
detachment rates 7(771) increasing with the mass m of the cluster where the particle is 
attached. Such dependence is certainly expected in heteroepitaxial growth if small lattice 
mismatches between the adlayer and the substrate are unfavourable to aggregation of new 
particles to existing islands. Indeed, the drastic consequences of this feature on island 
shapes was already illustrated in some systems, such as Cu islands on iVz(lOO), where a 
transition from compact to ramified shapes takes place as the coverage increases [20]. In 
one-dimensional systems, lattice mismatch is also expected to play a role, although at first 
approximation the shape effects are not present and the island mass is sufficient to determine 
the rate for detachment of a bordering atom. The association of a detachment rate with 



the island size may also be viewed as an alternative to describe t 



re effects of long-range 



22 



23 



24j, but which 



repulsive interactions present in a large number of real systems 21 
are usually very hard for computation. 

The models considered here are equivalent to zero-range processes and do have a fac- 
torisable steady state 25|, |26|. This feature means that the exact calculation of the steady 
state properties may be carried out with the use of an independent interval "approxima- 
tion" (IIA) to the master equation, which implicitly assumes a factorization of probabilities. 
The IIA, which has already proved to be an exceptional tool to investigate nonequilibrium 



statistical models 



19 



271 ]. here provides exact results for some forms of 7(771), which are 



confirmed by numerical simulation data. However, in many cases the latter approach is es- 
sential to calculate steady state properties, particularly when average cluster sizes are small 
and discretization effects play an important role. 



3 



We will show that both monomodal and monotonically decreasing cluster size distribu- 
tions may be obtained in these models depending on the particular form of detachment rate 
and the coverage. In the analysis of some forms of 7(777) which account for short-range 
attraction of neighboring atoms (via detachment rates much smaller than isolated particle 
diffusion coefficients), we will show that it is possible to exchange between those shapes 
by tuning the temperature or changing the coverage - high temperature and low coverage 
typically favoring the monotonically decreasing form. Although no quantitative comparison 
with real systems data will be shown here, we believe that the relative simplicity of our 
model and the range of qualitative behaviors obtained from it can motivate its use in par- 
ticular applications. From the theoretical point of view, this work opens the possibility of 
new applications of the widely studied zero-range processes. 

The rest of this work is organized as follows. In Sec. II we present the column picture in 
which the original problem (Fig. 1) is mapped, the general form of the master equation in 
the IIA and the method of solution in the steady state. In Sec. Ill we present the steady 
state cluster size distributions for selected forms of detachment rates and compare them with 
results of numerical simulations. We focus on the differences between limiting cases of rapidly 
and slowly increasing 7(771). In Sec. IV, we consider systems where formation of small islands 
is favoured by detachment rates significantly smaller than isolated atom (m = 1) diffusion, 
while 7(777.) is increasing. This mimics the competition between short-range attractive and 
long-range repulsive interactions, a case which may be especially relevant for surface science. 
In Sec. V we summarize our results and conclusions and discuss possible applications. 

II. PROCESSES, PICTURES, AND GENERAL FORMULATION 

As shown in Fig. la, the main processes in our problem are: (i) random walk of separated 
(single) particle, at a rate 7(1); (ii) detachment of particle from edge of cluster by particle 
stepping one unit away from it, at a rate 7(m), where m is the cluster mass. In addition, 
there is (iii) attachment of particle to edge of cluster, occurring immediately after the particle 
jumps to that position. 

The model defined in this so-called cluster picture can also be depicted in a column 
picture, as illustrated in Fig. lb. In the latter, the clusters are columns, and the group of n 
(> 1) vacant sites between two adjacent clusters has become a group of n — 1 (> 0) empty 
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columns between two filled columns. Thus, one of the vacancies of the cluster picture now 
acts as a column spacer, and the remaining ones have become columns with m = particles. 

This column picture corresponds to a zero-range process |25[]. Indeed the mapping just 
demonstrated achieves the equivalent of the inverse of a mapping from a zero-range process 



;o an exclusion process by other means 25|, |28(. In traffic flow models, such as that in Ref. 



29|], different rates for jumps to the left and to the right must be considered. 

The full analytic description of systems with such stochastic processes is provided by the 
Master equation, which is most easily written in the column picture. The description is 
simplified by the fact that the process conserves the total particle numbers N. Thus, using 
periodic boundary conditions and a total number of sites L (lattice length in the cluster 
picture), and denoting by N(m) the total number of clusters of size m (> 1), it follows that 
(i) iV = J2m=i m N(m), (ii) the number of spacers is J2m=i N(m), and (iii) N(m) equals the 
number of columns of size m, for m > 0. Hence, denoting by N(0) the number of columns 
of size zero, we have E^o^H = L — N = L(l — p) (the last step defining the density p 
in the original picture). Thus the total number of columns (including those of size zero) is 
constant, as is the density. The system configuration can be specified by the ordered set of 
numbers of particles in each of the columns in succession: (mi, . . .) = {m^}. 

The probability Pt{rrii} at time t of the configuration {m;} changes by in and out pro- 
cesses. For example, Fig. lc shows the process (mi, m 2 , m 3 , ...)—> (mi + 1, m 2 — 1, m 3 , . . .) 
(m.2 > 1) having rate 7(7712). Collecting the effects of all such processes in a time step 
t — > t + 1 gives the full Master equation 

L 

Pt+i{rrii} - P t {mi} = ^[7 (mj_i +l)P t (... mt-i + 1, m z - 1 . . .) 

+7 (m+i + 1) P% (• • • mi - 1, m l+1 + 1 . . .) 

-2j { mi ) P t {mi}}6 {mt) . (1) 

The theta function above (zero for m < 0, otherwise unity) is actually redundant as 
P (. . . m — 1 . . .) and 7(771) vanish for m < 0. 

As an ansatz, one can attempt to find a solution of the Master equation using for P t {mi} 
the factorised form Ilj=i Pti ( m i)- This turns out to give an approximate form (the IIA) for 
the time-dependent situation, but an exact result for the steady state (the time evolution 
may be particularly interesting in the case of rates decreasing with cluster mass, and is 
the subject of our current work). In the special case of homogeneous rates (7 independent 



of column position I) the steady state Master equation is solved exactly with the function 
P-a(m) — P{m), i. e. independent of t and I. Here, P(m) is the probability that an 
arbitrarily chosen column has occupancy m, i. e. N (m) /Em=o-^( m )i which corresponds 
to the probability that a cluster has mass m in the original problem (Fig. la). 

The reduced form of the steady state Master equation applying in this simplified (homo- 
geneous) situation can be obtained from Eq. ([1]) or directly as follows from the processes 
involved in the evolution of P(m). In an appropriately defined time step, the change of 
P(m) has positive and negative contributions (from in and out processes) 

P{m + l) 7 (m + l)0(m) + 6 mfi i(l)P(l) + T [P{m - l)9(m - 1) + S m ,iP(0)) (2) 

and 

- P(m) 1 (m)d(m - 1) - 5 m>l7 (l)P(l) - T [P(m)d(m) + 5 mfi P(0)] , (3) 
respectively, where 

oo 

r = 5:7(0^(0- (4) 

1=1 

The Master equation for the steady state now reduces to the condition (on P(m)) that 
the sum of contributions of in and out processes (Eqs. [2] and [3]) vanishes for each m. Defining 

A(m) = P(m) 7 (m) — TP(m — l),m > 1, (5) 

that gives the exact relations 

A(l) = 0, 

A(m + 1) - A(m) = , m > 1, (6) 
so that v4(m) vanishes for all m > 1. Thus P(m) satisfies 

r 

P( m ) = — — P(m- l),m > 1, (7) 
7 (m) 

yielding 

P(m)=P(0)n- m ,m>l. (8) 

Direct substitution of the resulting product form for P{m} into the steady state version of 
the original Master equation ([T]) verifies that this satisfies it exactly. This result is equivalent 



to the one 



261 ] previously given for zero-range processes (see e. g. Ref. [251] ). 
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Using the normalization condition 

oo 

£P(m) = l, (9) 

m=0 

it can be seen that Eq. ([8]) is consistent with the definition 
Using the definitions of Sec. II, we have 

(m) „ = —r ; — r- = > mr m = — : = , 10 

K UU J2™ =0 N(m) ^ { ] L(l-p) 1-p' 1 ] 

where (m) „ is the mean cluster size taking into account all columns, including those with 
zero mass. These equations provide a means to relate P(0) in Eq. (jSJ) to the density p, as 
well as the relationship between {m) all and p. 
Defining p = In T and 

m 

s(m) = In 7 (/), m > 1, 

s(0) = 0,m = 0, (11) 
and using Eqs. (|Hj) and we find 

oo 

P(m) = exp [m/i — s (m)]/ exp [i// — s (/)]. (12) 

z=o 

From Eq. (flQl) . /i can be determined in terms of p by 

= ^ mexp [m/x — s (m)]/ V] exp [m/i — s (m)]. (13) 

1 — P m=0 m=0 

The limits in the sums are those appropriate for an infinite system (N, L — > oo at fixed 
density p = N/L). The sums are constrained for finite L, N. 

For the infinite system there remains the question whether the sums converge or not. 
That depends on the form of the rates 7(771) at large m. In that region m can be treated 
as a continuous variable and sums become integrals. This continuum approach is also very 
useful for p near 1, where typical representative m's (such as {m) a u in Eq. [TO]) are large, so 
typically P(m) is appreciable at large m. Using such a continuum approach, it can be seen 
that for 7(777.) increasing with m at large m (making In 7 (m) positive and increasing), s(m) 
increases with m more rapidly than linearly, making the sums in Eqs. (fT2"|) and ffTB"]) 
converge, so the steady state solution for P(m) is physically acceptable: P(m) decreases 
with m at large m and is typically peaked. This is exactly the situation of interest in real 
systems where the increase of island size is unfavourable. 



In the continuum approach, the location m = mo of the peak in P(m) can be found using 



dm 



[fim — s (m) 



/i - hi7 (m ). 



(14) 



m=mo 



For mo 7' (mo) /7(mo) ^> 1 the peaking is strong, and P(m) can be approximated by 
exp \fx (m — m ) — / (m)], where / (m) = f™ In [7 {l)}dl. This has most of its 



■y(mo) 



2Try'(mo) 



1/2 



weight in the Gaussian form (P(m) oc exp — \ (7' (m ) /7 (mo)) (m — m ) 2 ), which applies 
near the peak. 

For 7(m) decreasing with m sufficiently fast at large m, the sums diverge and there is 
strictly no steady state in the infinite system. It can be shown that the criterion for no 
steady state in the infinite system is 7(m) decaying at large m more slowly than a constant 



times 1 + 



251 ] . In this case, the system will coarsen forever, but such situations [7 (m) 



decreasing with m] will not be considered here. 



III. STEADY STATE BEHAVIOR FOR SELECTED DETACHMENT RATES 

In this Section, analytic predictions and simulation results for steady state properties 
will be presented and compared. The emphasis will be on the empirically realistic "ther- 
modynamic" limit of very large systems, in which steady states are achievable. The steady 
state properties to be discussed here are the cluster size distributions and the average cluster 
sizes. However, it is important to mention that, from now on, cluster sizes are defined by 
averaging only over masses m > 1, in contrast with Eq. ( flOl) . which also took into account 
columns with zero mass. This new average will be denoted (m), and provides a more ap- 
propriate physical description of the system in the original cluster picture. Indeed, except 
where explicitly indicated, we will refer to that original picture in the following. 

The simulations were typically performed in lattices of sizes L = 8192, with several den- 
sities and different forms of 7(m). Due to the small cluster sizes imposed by the system 
dynamics, in all cases finite-size effects are negligible (this was checked by comparison of 
results in different lattice sizes). The generation of a sequence of configurations begins 
with the deposition of a random layer of density p. Subsequently, the dynamics with diffu- 
sion, attachment and detachment processes is allowed, and the evolution of the cluster size 
distribution is monitored. Typically, it is assumed that the steady state is attained if no 
appreciable change (e. g. 1% in the peak) is found in the distribution at the last time decade 



8 



of the simulation, after averaging over at least 100 initial configurations (i. e. 100 different 
sequences). Under these conditions, the representative configurations are attained, on the 
average, after 100 — 1000 detachments of particles aggregated to small clusters. Anyway, it 
is important to note that in our models steady state properties are independent of the initial 
system configuration, thus the same results would be obtained if diffusion and detachment 
processes were competing with deposition during the production of the first configuration 
with the desired density p (however, there are special models where steady state properties 
depend on the initial configuration - see e. g. Ref. |32j). 

First we consider the case (i) of constant detachment rate, 7(m) = a for m > 2 and 
isolated particle hopping rate 7(1) = b. The analysis proceeding from Eqs. (1131) and (IT2|) 
gives a size distribution which is exponentially decreasing in m for m > 2. For the special 
case b = a, using s(m) = mho. a, we have P(m) = (1 — p) p m for all m. Such results were 
formerly predicted in Ref. jl9j] and numerically confirmed in Ref. 30]. These distributions 
are equivalent to those obtained in equilibrium (reversible) models with nearest neighbor 



interactions between the adatoms 



H3|. 

Now we consider the power law case (ii), 7(771) = m k for m > 1, in which the reduction of 
the analytic result (fT2l) for P(m) is less straightforward than in the previous case. For this 
reason, most results are obtained from simulation. In Fig. 2, we show the scaled distributions 
for k = 1 and densities p = 0.5, p = 0.75 and p = 0.95 (in this and subsequent plots, dashed 
curves are guides to the eye). For small and medium densities, isolated particles (m = 1) 
are predominant and the distribution is rapidly decreasing, so that it seems to decay faster 
than a simple exponential for p < 0.5. As the coverage increases, it crosses over to a peaked 
(monomodal) distribution, which becomes very sharp for p close to 1. This crossover is 
directly related to the increase in the average size (m) « pj (1 — p), from a value near unity 
for small p (where the decrease is monotonic) to large values for large p (where the peak is 
close to (m)). 

The analytical distribution of case (ii) becomes simple for p ~ 1, since there the typical 
masses are large and continuum approximations can be used. That results in P(m) oc 
exp [(p + k)m — fcmlnm] (m large). For k — 1, this is the large m approximation to a 
Poisson distribution. The agreement with simulation results is illustrated in the inset of 
Fig. 2 for k — 1 and p = 0.95, where log [P (m)} is shown to be a function of m In (m) — Cm, 
with a fitting constant C = 3.92. 
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We next consider another simple detachment rate function, 7(771) = (^pj) f° r rn > 1, 
hereafter called case (iii). The main difference from case (ii) is the fact that 7 does not diverge 
asm-* 00. Using Eqs. ({TBI and ( fl2l) . we obtain exactly P{m) oc (m + l) b exp {—(3m), with 
(3 constant. For low integer values of b, the constant (3 and the normalization constant can 
easily be obtained analytically in terms of the density. For the simplest case h = 1, we obtain 
j3 = In i^r) an d (771) = (1^-^4-3^ ■ We remark here that the case b = —2 is marginal: from 
the discussion at the end of Sec. II (see also Ref. |25j), the infinite system achieves a steady 
state if b > — 2 but not if b < -2. 

In Fig. 3 we show the scaled cluster size distributions for case (iii) obtained from sim- 
ulation, with 6=1 and coverages p = 0.5 and p = 0.9. Both show excellent agreement 
with the analytical results (hereafter represented by solid curves in the plots). Again we 
observe that the distribution is monotonically decreasing for small coverages, where (m) is 
close to 1, while for larger coverages there appears a peak. However, there is an important 
difference from case (ii) here: instead of having a sharp peak close to (m), the distribution 
for large coverages has a very fat left tail and the most probable island size may be smaller 
than (m)/2 (see data for p = 0.9 in Fig. 3). This difference is certainly a consequence of 
the slower increase of the detachment rate with the cluster mass. 

We also analyzed other forms of detachment rates, namely simple exponential, logarithmic 
7(m), and the particular logarithmic case 7(771) = c[ln (m + i)][ 1+ln ( m+1 )l ; f or which a closed 
analytical form of P(m) can be obtained. In all cases, the analytically predicted distributions 
agree very well with simulation data. With the logarithmic forms, the qualitative behavior 
is intermediate between cases (ii) and (iii) above. 

The peaked cluster size distributions shown above were always obtained for large densities. 
At first sight, this condition seems to be unattainable in experiments on sub monolayers 
on vicinal surfaces, since the one-dimensional character of the adlayer (e. g. chain-like 
structures along the step edges) is lost at high coverages. However, in Sec. IV, we will 
show that under certain realistic conditions on 7(m), it is also possible to observe peaked 
distributions for densities lower than 0.5. 

Moreover, it is important to stress that the density p of particles along the steps can 
be much larger than the nominal coverage 9 of the surface because the latter is an adatom 
density per substrate area, while p is the filling of one- dimensional rows. As an example, 
we refer to the STM image of Ag deposited on Pt(997) in Fig. lb of Ref. [6|]), where the 
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coverage 6 = 0.04 monolayers is much smaller than the effective filling p fa 0.3 of the step 
rows. 



IV. SYSTEMS WITH COMPETING INTERACTIONS AND TEMPERATURE 
EFFECTS 

Here we extend the study of Sec. Ill to systems with the mass-dependent detach- 
ment rates all far less than the isolated atom diffusion rate 7(1), by an overall factor 
e ~ exp (—E/kBT) <C 1. Here, E may be interpreted as a binding energy between neighbor- 
ing adatoms, and T is the temperature (see e. g. the discussion on the diffusivities for initial 



state interactions in Ref. |3l|], where only short-range interactions were considered). For 
low temperatures, that factor significantly reduces the mobility of aggregated atoms when 
compared to the isolated ones. On the other hand, the mass dependence of 7 accounts for 
the effects of the interactions with the substrate, which works against the increase of cluster 
size (e. g. effects of lattice mismatch or substrate-mediated repulsive interactions). We will 
typically work with 10 _1 < e < 10~ 3 for each form of 7(771), in order to show the possible 
effects of the temperature in the island size distributions. Simulation work here will focus 
on low coverages, typically below or at half filling of the one-dimensional rows. 

First we consider the generalization of case (ii), where 7(771) = c(m)m h with c(l) = 1, 
and c(m) = e < 1 for m > 2. This means that detachment rates from small clusters (up to 
masses ~ e -1 ) are smaller than the isolated atom diffusion rate, but larger clusters are very 
unstable. 

For k = 1, no simple closed form for the cluster size distribution can be obtained. In Figs. 
4a and 4b we show the simulation results with densities p = 0.1 and p = 0.5, respectively. In 
both plots we consider e = 10 -2 and e = 10~ 3 . For low density (Fig. 4a), the distributions are 
monotonically decreasing up to e ~ 10~ 2 , but a peak at small m appears at sufficiently low 
temperatures (i. e. very small e). For medium density, Fig. 4b shows that the temperature 
does not need to be so small for the onset of a peaked distribution: the peak is present for 
e = 10 -2 and it is well defined at e = 10~ 3 . These results must be compared with those in 
Fig. 2 for k = 1 and p = 0.5, but e = 1, where the distribution is monotonically and rapidly 
decreasing. 

Analytical results can be obtained for the generalized case (ii) only for k < 1 and e <C 1, 
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where average cluster sizes are large and a continuum approximation of the cluster size distri- 
bution is possible. In these conditions, we have P (to) oc In B/(\nB)^ } , with B = 
( J^ ) (i-k)/2 e (i-p) an( i C — ^p- This is a monomodal distribution with (to) « | In _B/(ln_B) < ' , 
which confirms the general trend that this shape is favoured by large densities and small 
temperature (small e). 

Unfortunately, the above formula for P{m) is accurate only for very small e and very 
small k, otherwise the continuum approximation fails due to the discreteness of the typical 
cluster sizes and the large statistical weight of isolated particles. As an example, we compare 
in Fig. 5a and 5b the analytical and numerical cluster size distributions for k = 1/2 and 
p = 0.5, with e = 10~ 2 and e = 10~ 3 , respectively. Even for e = 10~ 3 , where (to) ~ 7, we 
observe deviations of the analytical approximation from the numerical data. 

The results in Figs. 4 and 5, as well as the analytical approximation for k < 1 and 
e< 1, show another important feature: the peaks of the distributions are very close to (m), 
similarly to the case (ii) with e = 1 studied in Sec. III. We recall that these are cases of 
rapidly increasing 7(m), in which the formation of large clusters is highly unfavourable. 

Now we consider the generalization of case (iii), where 7(771) = c (to) (j^ij , with c(l) = 
1, and c(m) = e for to > 2. The resulting probability distribution P(m) is of the same 
form as the one in Sec. Ill (i. e. for the special case e = 1), however the small e gives 
much lower weight to P(m), to > 1, than P(0), and the peak of the distribution is, for most 
densities, pushed out to much larger masses. For general 6 it is possible to show that for 

small enough e (very much less than both p and (1 — p) b+1 ) the mean cluster size has the 

r (b+i) b + lll/{b+2) 
form (771) - — - — - 



_e(l-p) b\ 

In Figs. 6a-d we show the simulation results for 6 = 1, with densities p = 0.1 and p = 0.5 
and detachment factors e = 10~ 2 and e = 10~ 3 . For this value of 6, the full distribution can 

1 /3 

be obtained analytically as P(m) oc (m + 1) e _/3m , (3 = ( 2e ^~^ ) • It compares well with 
simulation data for intermediate densities and small values of e, as shown in Fig. 6b. 

Lowering the temperature also provides peaked distributions in this case, and the position 
of the peak is shifted to larger masses as the density increases. The effect of temperature is 
important only in the right tail of the distribution, which typically has a small number of 
data points, corresponding to small islands. Similarly to what was observed in Sec. Ill, the 
peaks are located close to (m)/2 and the left tail is very fat. This is a signature of a weak 
size-dependence of the detachment rates. However, the main difference from the results in 
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Sec. Ill is again the possibility of finding the transition from a monotonically decreasing 
distribution to a monomodal one by tuning the temperature, with a coverage not too large. 

For the general case of rate functions 7(m) increasing with m (at large m), the mean 
cluster size (m) can always be made large at low densities by suppressing detachment rates 
by a suitably small factor e, compared to the diffusion rate 7(1). Considering 7(m) = tv(m) 
for m > 2 and 7(1) = ^(1), the relationships 

{m)~F(X) , X = -^-^ (15) 

of average cluster size to density and Arrhenius rate parameter e which result at small e for 
a variety of detachment rate functions (at large m) are: 

u(m) = m k F(X) = ^mA, 

rC 

Km) = (lnmr ^ F(x) = _j^_, 



. In .V 



v{m) = [me [lnm > \ F(X) = 
v(m) = e am " F(X) = 



a In In X 

an 



l/(n+l) 



(16) 



It can be seen that the effect of e (in reducing the density required to get a large mean 
cluster mass) diminishes as the rate of increase of v with m increases. More importantly, 
these results show that density and temperature control the island size distribution through 
the scaling variable X, so that monomodal (monotonically decreasing) forms are found for 
large (small) X. 



V. DISCUSSION AND CONCLUSION 



We studied one-dimensional models of island formation by diffusion, attachment and 
detachment of single particles, considering detachment rates 7 increasing with the island 
mass m. This type of model may be of experimental interest because it can be used to 
simplify the description of more complex interactions that prevent the formation of large 
atomic islands. Indeed, these models are equivalent to zero-range processes in which the 
system attains a steady state. An independent interval "approximation" to the Master 
equation was proposed in order to calculate the cluster size distribution which actually 
provides an exact description in the steady state. The tools necessary for derivation of 

13 



explicit distributions for any particular form of 7(771) were provided and results for a variety 
of cases were presented and compared with numerical simulation data. 

The representative rate functions analyzed above include some which arise from associat- 
ing (Arrhenius) detachment rates with potentials U(m) for particles at the end of a cluster 
of size m. U(m) is then a sum, from I = 1 to m — 1, of pair potentials V(l) for separation 
I. The forms V vanishing, V Coulomb-like and inverse square then give (exactly in the first 
two cases, and for large m in the third one) the detachment rate functions of cases (i), (ii), 
(iii), respectively. In addition, we also mimicked the presence of short-range attraction in 
the models where the detachment rates were smaller than isolated particle hopping rates by 
an overall factor e < 1. 

One of the important conclusions of this work is the possibility of changing the shape 
of island size distributions, from monomodal to monotonically decreasing ones, by tuning 
the temperature or changing the coverage. For different forms of detachment rates (rapidly 
or slowly decreasing with m), it was shown that high temperatures and low coverages are 
favourable for the onset of the monotonically decreasing distributions, their combined effect 
being described by the scaling variable X = ^7^3 ■ This transition is not usual in non- 
equilibrium modeling of submonolayer growth, particularly when irreversible attachment to 
islands (i. e. a critical island size) is assumed. However, it is an important finding because 
both shapes of distribution were already observed (isolated) in island growth near step edges 
of vicinal surfaces [5|, |6| . 

On the other hand, rapidly or slowly increasing 7(771) have different effects on the position 
of the peak of the island size distributions and on its right tail. Rapidly increasing 7(771) 
[e. g. case (ii) above] lead to peaks very close to (m) and rapidly decreasing tails, since 
formation of large clusters is very difficult. However, the slowly decreasing 7(7/7) forms [e. 
g. case (iii), where 7(777) does not diverge when m — > 00] provide peaks close to (m)/2 and 
fat right tails. In real systems where island size distributions are measured, these results 
may give clues on how intense are the mechanisms that work against the formation of large 
clusters. 

Despite the fact that no application to a particular real system was proposed here, we 
believe that the framework developed in this paper may be useful for such applications. Long 
range interactions frequently play a role in submonolayer growth but introduce difficulties 
to both analytical (scaling) and numerical calculations, even when they are limited to one 
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dimension (see. e. g. 33]). Consequently, the inclusion of simpler mechanisms in a model 
system may be useful, such as the association of detachment rates to the full cluster size 
suggested here. The fact that our model corresponds to a zero-range process allows for 
much simpler analytical calculations of steady state properties, in contrast to other (not less 
important) approaches, such as the introduction of energy barriers for particle attachment 



to clusters 



23 



34j . Finally, it is also interesting to recall that gas adsorption in carbon 



nanotubes may be viewed as a one-dimensional clustering problem, thus it is another field 
where this type of non-equilibrium model may find application. Indeed, simple statistical 
equilibrium models of interacting particles in finite lattices were already proposed for those 



systems 



35]. 
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FIG. 1: (a) Illustration of the diffusion (m = 1) and detachment (m > 1) processes of the model, 
with the associated rates 7(171). (b) Examples of detachment processes (1,2) and hopping processes 
(3,4) of filled particles, in the original cluster picture and in the corresponding column picture. 
Dashed lines show the correspondence between clusters in the two pictures, (c) Example of leftward 
movement of an aggregated particle. 
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12 3 

m/ <m> 

FIG. 2: (Color online) Scaled cluster mass distribution for detachment rate 7(777) = m (case (ii) 
with k = 1), with densities p = 0.5 (squares), p = 0.75 (triangles) and p = 0.95 (crosses). Dashed 
curves are drawn to guide the eye. Inset: the rescaled data for 7(771) = m and p = 0.95, with a 
fitting constant C = 3.92, is consistent with a Poisson distribution (the solid line is a linear fit of 
the data). 
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12 3 

m/ <m> 

FIG. 3: (Color online) Scaled cluster mass distribution obtained from simulation for detachment 
rate 7(m) = m/(m + 1) [case (iii)] with densities p = 0.5 (squares) and p = 0.9 (triangles). Solid 
curves are the corresponding analytical results. 
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FIG. 4: (Color online) Scaled cluster mass distributions obtained from simulation for detachment 
rate 7(771) = c(m)m, with c(l) = 1 and c(m) = e for m > 2, and densities (a) p = 0.1 and (b) 
p = 0.5. In both panels, e = 10 -2 (triangles) and e = 10~ 3 (squares). Dashed curves are drawn to 
guide the eye. 
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12 3 
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FIG. 5: (Color online) Scaled cluster mass distributions for detachment rate 7(771) = c (m)m 1 / 2 , 
with c(l) = 1 and c(m) = e for m > 2, density p = 0.5 and (a) e = 10 -2 and (b) e = 10~ 3 . 
Symbols indicate simulation data and solid curves show the corresponding analytical results in the 
continuum approximation. 



22 



1 

0.8 
0.6 

4 

3 0.2 
Oh 

0.8 

0.6 

0.4 

0.2 





A 

s 

V 



1 1 1 1 1 1 1 1 1 1 1 1 1 

r'X (a) : 


i i i i 1 i i i i 1 i i i 

(b) ~ 


1 1 1 1 | 1 1 1 1 | 1 1 1 

7 r\ (c) : 

i i i i 1 i i i i 1 i i i 


1 1 1 1 1 1 1 1 1 1 1 1 1 
(d) " 

1 1 1 1 1 1 1 1 1 1 1 1 1 



12 12 

m/ <m> 



FIG. 6: (Color online) Scaled cluster mass distributions for detachment rate 7(771) = 
c(m)m/(m + 1), with c(l) = 1 and c(m) = e for m > 2, and: (a) p = 0.1, e = 10~ 2 ; (b) 
p = 0.1, e = 10~ 3 ; (c) p = 0.5, e = 10~ 2 ; (d) p = 0.5, e = 10~ 3 . Symbols indicate simulation data 
and solid curves show the corresponding analytical results. 
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